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Abstract 

Motivated by a number of recent experimental studies we have revisited the problem of the mi- 
croscopic calculation of the quasiparticle self-energy and many-body effective mass enhancement in 
a two-dimensional electron liquid. Our systematic study is based on the many-body local fields the- 
ory and takes advantage of the results of the most recent Diffusion Monte Carlo calculations of the 
static charge and spin response of the electron liquid. We report extensive calculations of both the 
real and imaginary parts of the quasiparticle self-energy. We also present results for the many-body 
effective mass enhancement and the renormalization constant over an extensive range of electron 
density. In this respect we critically examine the relative merits of the on-shell approximation 
versus the self-consistent solution of the Dyson equation. We show that in the strongly-correlated 
regime a solution of the Dyson equation proves necessary in order to obtain a well behaved effec- 
tive mass. The inclusion of both charge- and spin-density fluctuations beyond the Random Phase 
Approximation is indeed crucial to get reasonable agreement with recent measurements. 
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I. INTRODUCTION 

An interacting electron gas on a uniform neutralizing background (EG) is used as the 
reference system in most realistic calculations of electronic structure in condensed-matter 
physics yj. At zero temperature there are only two relevant parameters for a disorder- free 
EG in the absence of quantizing magnetic fields and spin-orbital coupling: (i) the usual 
Wigner-Seitz density parameter r^ = (7^20^5) ""^^^j '^b = ^^^/(me^) being the Bohr radius 
in the medium of interest with R and m appropriate dielectric constant and bare band mass 
respectively; and (ii) the degree of spin polarization ( = \n^ —n^l/n^u- Here Uo- is the average 
density of particles with spin a =t, i and naD = ""-t + '^i is the total average density. 

Understanding the many-body aspects of this model has attracted continued interest for 
„.any decade. M The EG, unlike .y=te„. of cla..cal pa.t.cte, behave, hke an ideal 
paramagnetic gas at high density (r^ <^ 1) and like a solid at low density |^ (r^ ^ 1). In 
the intermediate density regime, which is relevant in three dimensions (3D) to conduction 
electrons in simple metals and in two dimensions {2D) to electrons in an inversion layer 
of a Si metal-oxide-semiconductor field-effect transistor (MOSFET) or in an AlGaAs/GaAs 
quantum well, perturbative techniques are not effective owing to the lack of a small expansion 
parameter. One has to take recourse to approximate semi-analytical methods, a number 
of which have been reviewed in Ref. M, or to Quantum Monte Carlo (QMC) simulation 



methods 
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Among the methods designed to deal with the intermediate density regime, of particular 
interest for its physical appeal and elegance is Landau's phenomenological theory |ly| dealing 
with low-lying excitations in a Fermi-liquid. Landau called such single-particle excitations 
quasiparticles (QP's) and postulated a one-to-one correspondence between them and the 
excited states of a non-interacting Fermi gas. He wrote the excitation energy of the Fermi- 
liquid in terms of the energies of the QP's and of their effective interactions. The QP-QP 
interaction function can in turn be used to obtain various physical properties of the system 
and can be parametrized in terms of experimentally measurable data. 

Quinn and Ferrell [17[ provided a framework for the microscopic evaluation of the QP- 
QP interactions in the EG by means of the Random Phase Approximation (RPA). Next, 
Rice Il8j incorporated the vertex corrections in the RPA form of the electron self-energy 



by including the Hubbard [19| many-body local-field. Some problems in Rice's theory were 



subsequently resolved by Ting, Lee, and Quinn [20| in their theory of the quasi-2D EG. All 
these approaches considered only the charge-density fluctuations while neglecting the effect 
of the spin-density fluctuations. 

A more detailed analysis that accounts for the vertex corrections associated with both 



types of fluctuations was carried out for an unpolarized EG in Refs. 



21 



m 



23, where 



Kukkonen-Overhauser-like [2J] effective interactions were obtained by different approaches. 
In particular, Yarlagadda and Giuhani 



I adopted a physically transparent approach 
termed renormalized Hamiltonian approach (RHA), which will be extensively discussed in 
this paper. A few electrons from the EG are selected and called "test electrons", while 
the remaining EG is treated as a dielectric screening medium. As the test electrons move 
through this medium, they produce fluctuations in the density of spin-up and spin-down 
electrons, which provide virtual clothing and also screen their interactions. Thus, the di- 
electric mimics the true physical processes in an average way. Of course, the test electrons 
and the electrons of the medium are physically indistinguishable, and this must be taken 
into account when exchange effects are considered. At this point, after averaging over the 
coordinates of the screening medium, an effective renormalized Hamiltonian containing only 
the degrees of freedom of the clothed test electrons (or QP's) can be derived, under the 
assumption that the coupling with the medium occurs only via its charge- and spin-density 
fluctuations. The basic idea underlying the RHA was earlier developed in a beautiful paper 
by Hamann and Overhauser 12511 within the RPA. Calculations based on these theories have 

Tin ^ 

been carried out for both 3D |22, |2^ and 2D [27] systems. 

In a parallel theoretical development Ng and Singwi [2^, starting from a Ward identity 
and performing a local approximation on the irreducible particle-hole interaction, obtained 
an expression for the self-energy in terms of the many-body local field factors associated 
with charge- and spin-density fluctuations. Equivalent results were later obtained by Yarla- 
gadda and Giuliani (29l . |2i| by means of the RHA. These authors also took into account an 
infinitesimal degree of spin polarization which allowed them to carry out calculations of the 
Landau Fermi-liquid parameters. 

The previous theoretical work suffers from two major shortcomings, (i) All available 
earlier calculations have adopted a static and oversimplified Hubbard-like model for the local- 
field factors, which do not have the appropriate behavior at both intemediate and large wave 
number q. We have corrected for these discrepancies (while still neglecting the frequency 



dependence of the local fields) using recent parametrizations 



for the static local fields 



of a 2D EG. (ii) The RHA briefiy introduced above is an "on-shell" theory (see Sect. IIIB|) 
and, as we will show in this work, predicts a divergence of the effective mass with decreasing 
electron density. In the spirit of the work of Santoro-Giuliani J27l] and of Ng-Singwi J28J we 
have kept in this work the full frequency dependence of the self-energy and carried out a 
self-consistent solution of the Dyson equation to find the proper QP excitation energy and 
QP properties. Comparing with Ref. uA we have released the plasmon+paramagnon-pole 
approximation to the charge-charge and spin-spin response functions. 

From the experimental point of view, as already remarked, electrons in a semiconductor 
inversion layer or in a quantum well can be modelled by a quasi-2Z} EG [32j. Quantum 
Shubnikov-de Haas (SdH) oscillations of the magnetoresistance 33| provide a powerful tool 



for measurin 
past years J34 



ermi-liquid parameters of a quasi-2Z} EG. Measurements performed over the 
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37[ have shown sizeable renormalizations of the QP effective mass and 
effective Lande g-iactoi. These experiments have been performed in a relatively high-density 
regime, i. e. for r^ < 2 say. The density dependence of the effective mass m* was obtained 
by Smith and Stiles 134 from a study of SdH oscillations in Si inversion layers. To obtain 
the same information Abstreiter et al. |35|] used instead cyclotron resonance measurements. 
Fang and Stiles |3fl] and Neugebauer et al. J37[ performed a series of SdH experiments on 
Si inversion layers and obtained the dependence of the modified Lande factor g* on carrier 
density. The product of g* and m*, which is proportional to the spin susceptibility xs, can 
be determined from the SdH oscillations in a tilted magnetic field as suggested in Ref. I3(tI . 

The issue of the apparent metal-insulator transition [2q] (MIT) in low-density 2D electron 
systems has prompted intense experimental studies on quasiparticle properties [39|, |40|, |4l|, 
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46l |47[ in the intermediate-to-strong coupling regime, r^ > 2 say. Many 



39'] have shown that the resistance of a Si-MOSFET is increased dramatically by 



increasing the value of an in-plane magnetic field, and saturates at a characteristic value 
of several Tesla. Performing low-field SdH measurements on Si-MOSFET's, Okamoto et 
al. J4i| have shown that the saturation value is the r nag netic field that is necessary to fully 



42| | of the in-plane magnetoresistance 



polarize the electron spins. An interpretation 4l|, 

in Si inversion layers suggested a ferromagnetic instability at or very close to the critical 

density for the 2D MIT driven by a divergence in the effective mass. Direct measurements 



of m* in high-mobility Si-MOSFET's over a wide range of carrier density, using a novel 



technique based on the beating pattern of SdH oscillations in crossed magnetic fields, have 
been reported by Pudalov et al. 43|. These authors measured m* and xs in the vicinity of 
the 2D MIT, but found no evidence for a divergent behavior. Only a moderate enhancement 
of m* by a factor of ~ 2—2.5 over the band mass was observed near the critical density for the 
2D MIT. Two groups have also reported anomalous density dependences of the Lande factor 
in n-doped J4^ (2 ^ t^ < 7) and p-doped 45| (r^ > 17) GaAs/AlGaAs heteroj unctions that 
are in disagreement with results in Si-MOSFET's. The dependence of the spin susceptibility 
on the degree of spin polarization of the sample can account for this anomalous behavior as 
pointed out by Zhu et al. |4a| , who studied a 2D EG of exceedingly high quality. 

To complete the cornucopia of recent experimental findings on QP properties, it is worth 
mentioning that Vakili et al. |47j have reported measurements of m* and xs iii a dilute 
2D EG confined to a narrow AlAs quantum well (only 45 A wide). The electron system 



investigated in Ref. |47| is quite interesting because the electrons occupy an out-of-plane 
conduction-band valley, rendering the system similar to 2D electrons in Si-MOSFET's but 
with only one valley occupied. Quite surprisingly, the results of Vakili et al. [42] for xs 
are in good agreement with the QMC results of Attaccalite et al. |l^ even though this 
simulation has been carried out for a strictly disorder- free EG. This might indicate that xs 
is not strongly dependent on disorder. On the other hand, there is a significant spread in 
the experimental results of Ref. |47| for m* , which turns out to be both sample and cool- 



down dependent. Difficulties associated with the SdH data analysis have been pointed out 



in Ref. |42| as one of the possible causes for this spread. 

At this point it is probably worth commenting that indeed there could be in principle 
subtle issues associated with the analysis of the SdH traces in 2D systems. In fact the ampli- 
tude of the SdH oscillations is usually fitted to the 3D Lifshitz-Kosevich (LK) formula 48| 
upon a trivial change in the single-particle spectrum. The fit is based on an impurity 
scattering Dingle temperature To and an "effective" mass. In recent years a number of 
caveats concerning the applicability of such a procedure to 2D strongly-interacting systems 
have appeared |42|. In particular, Martin et al. pj] have shown that the interplay between 
electron-electron interactions and electron-impurity scattering leads in 2D to an effective 
temperature-dependent Dingle temperature with a leading low-temperature behavior of the 
type Td{T) oc TlnT. The need for the introduction of a temperature-dependent Dingle 
parameter in strongly-coupled Si-MOSFET's has been emphasized in the above-mentioned 



Ref. |43[ where a linear Td{T) was used to fit the longitudinal magnetoresistance data. 
Quantitative differences on the resultant effective mass are found using such type of proce- 
dure: roughly speaking, the tendency is to get substantially lower values for m* than those 
obtained using the same Dingle parameter for all temperatures. 

For a quantitative comparison between suitable theories that take into account quasi-2D 
effects (such as finite width of the electron wavefunctions in the confinement direction and 
valley degeneracies) and the experimental results J3J, |36[ for Si-MOSFET's in the weak- 
coupling regime r^ < 2 we refer the reader to the work of Yarlagadda and Giuliani 30| and 
references therein. In this work we will trv and carry out a comparison between the theory 
and the experimental data of Zhu et al. J51| for strongly-interacting electrons (2 < r^ < 6), 
occupying a single valley in an exceptionally clean GaAs/AlGaAs quantum well. 

The contents of the paper are described briefiy as follows. In Sect. |n]we present in great 
detail the theoretical background. We proceed in Sect. lIIII to discuss the inputs we have used 
for our numerical calculations, while in Sect. II VI we present our main results for the real and 
imaginary part of the quasiparticle self-energy, the many-body enhancement of the effective 
mass, and the renormalization constant. Finally, in Sect. |3 we compare our theory with 
the experimental results of Zhu et al. [51j and report some conclusions. In order to make 
the paper fully self-contained we have included two Appendices which contain a number of 
helpful details on how we have in practice calculated the QP self-energy. 

II. THEORY OF THE QUASIPARTICLE SELF-ENERGY 

The aim of this Section is to give a sound theoretical justification to the expression for the 
retarded QP self-energy Si.ot(k, i^) of a 2D paramagnetic EG that will be used throughout 
this paper and that we have summarized in Eqs. Q and ((Tj) below. 

To fix the notation we start by introducing the Hamiltonian for a 2D EG confined to an 
area *i? , 

^EG = Y^ £kci,<xCk,a + l^Y^'nY^ Y^ 4i+q,<xi42-q,<X2Ck2,a2Cki,ai ■ (l) 

k,cr q^O ki,(Ti k2,o-2 

Here Cj[ ^ and Ck.cr are fermionic creation and annihilation operators which satisfy canonical 
anticommutation relations, e\^ = h^\^ /{2m) is the single-particle energy, and v^^ = 27re^/g 
is the 2D Fourier transform of the bare Coulomb interaction e^/r. For later purposes 



we introduce the Fermi wave number kp = (27^20)^'^^ = V2/{rsaB), the Fermi energy 
ep = h'^kp/{2m) and the quantity C,k = s^ — ep- 

The retarded QP self-energy Sret(k, u;) is written as the sum of two terms, 

Eret(k, u) = Ssx(k, u) + ScH(k, uj) (2) 

where the first term is called "screened-exchange" (SX) and the second term is called 
"Coulomb-hole" (CH). The frequency u is measured from ep/h. 
The SX contribution is given by 

Ssx(k, ^) = - / yf^ -, V^M 0(-^kW^) • (3) 

J {2-kY e(q, w - ^k+q/^) 

Here Q{x) is the step function and eia^uo) is a screening function originating from effective 
Kukkonen-Overhauser interactions [2J|, 

l + v^ {l-G+{(i,u)f xc(q,^) + 3t;qG^(q,^)Xs(q,^)- (4) 



In Eq. (0)) the charge-charge and spin-spin response functions xc{^.i^) ^"^^ XsIQ;^) ^"^^ 
determined by the spin-symmetric and spin-antisymmetric local-field factors G'+(q, cu) and 
G_(q,^), 



Xcq,w =^ r, pTT T^ — 7 V (5) 

l-i;q[l-G'+(q,^)]Xo(q,^) 

1 + i;qG'_ q, a; xo q, w 



and 



Xo(q, 1^) being the Stern response function of a noninteracting 2D EG J5^. In the param- 
agnetic electron liquid G±(q, cj) = [(^^^(q, cu) ± G'|j(q, co')]/2, where (^^^'(q, ^) are the spin- 
resolved local fields. Note that Ssx(k, cu) is just an ordinary exchange-like self-energy built 
from the Kukkonen-Overhauser effective interactions instead of bare Coulomb interactions, 
which would lead to the frequency-independent Hartree-Fock self-energy first calculated for 



the 2D EG by Stern [53 1. 



The CH contribution to the retarded self-energy is given by 

EcH k, w = - / — --^ i;q / — p— — T, (7) 



where (5 is a positive infinitesimal. The real and imaginary part of the retarded self-energy 
are readily obtained from Eqs. and ((7j) with the result 

3fJeSret(k, ^)= - j -0-2 ^q3f^e[£-i(q, uj - ^k+q/^)] e(-a+q/^) (8) 



(27r)2 'I io TT cu-aW^-^' 

and 

3mSrct(k,a;) = / — -^ i;qS5m[£"^(q,^ - ^k+q/^)] [©(^ " ^k+q/^) - ©(-^k+q/^)] • (9) 

Once the QP self-energy is known, the QP excitation energy (5£^Qp(k), which is the QP 
energy measured from the chemical potential /i of the interacting EG, can be calculated by 
solving self-consistently the Dyson equation 

5^Qp(k) = ^k + ^eS^,(k,^)|_,^^^(^^/^ , (10) 

where 3fJeEj^^(k, a;) = 3fJeSi.ot(k, cj) — Srctl^F, 0). For later purposes we introduce at this point 
the so-called on-shell approximation (OSA). This amounts to approximating the QP excita- 
tion energy by calculating 3fJeS^j(k, oj) in Eq. (fTn|) at the frequency u = ^k/h corresponding 
to the single-particle energy, that is 

<5£Qp(k) ^a+ 3?eSfet(k,^)L^^^/, . (11) 

In Sects. Ill Al and III Bl we will give a formal justification of Eqs. and ((2j) using two 
completely different methods: diagrammatic perturbation theory fSect. ITT^ and renormal- 
ized Hamiltonian approach (Sect. I11B|) . In the following we shall soon drop, however, the 
frequency dependence of the local-field factors. Recent studies J5J| have evaluated it in the 
long- wavelength limit q ^ 0, but the knowledge of the full dependence on wave number is 
necessary for the type of calculations that we are interested in. 

A. Diagrammatic perturbation theory 

In this Section we use a diagrammatic approach that was first developed by Ng and 
Singwi 



2811 ■ building on earlier ideas by Vignale and Singwi ^]. The starting point is the 

-7^^/k..,k'..v'(0)<5a'(k', u;') , (12) 



identity [5 



6G^,(p') 




5la(P) = 



FIG. 1: Diagrammatic representation of the identity dJ). Here and in the fohowing we use the 
four-momentum variable p as a shorthand for (k, lj). 

where 5So-(k, cj) and 5Go-(k, cj) are infinitesimal changes in the self-energy and the Green's 
function, and Ikuja,ii'uj'a'{^) is the irreducible electron-hole interaction at zero momentum and 
energy transfer. This identity is graphically represented in Fig. 1. The defining feature of 
the irreducible electron-hole scattering block / is that it includes only diagrams that cannot 
be divided into two parts by cutting a single electron-hole pair propagator carrying zero 
energy and momentum. 

The differential relation ()12p cannot be integrated as it stands because / is a complicated 
functional of G. The idea of Ng and Singwi was to use an approximate form of /kwa,k'w'cr'(0) 
that does not depend on G. The "local approximation" introduced by Vignale and Singwi 
in their study of the effective electron-electron interaction ^] is useful for this purpose since 
it yields by physical arguments an expression of the form 

/k..,kw(0) ~r;5(k-k',a;-cu') , (13) 

where V^^, is just a function of the momentum and energy transfers in the electron-hole 
channel. Thus the main characteristic of the Ng-Singwi approach is that the key approxi- 
mation in Eq. M'J\) is made on the irreducible electron-hole interaction rather than on the 
self-energy itself. With this approximation we can integrate Eq. ()12|1 and obtain, up to an 
integration constant, the result 

^-^VfA\.-\^^u-J)GA^^u') . (14) 

With the replacements k— k' = q and uj—uj' = Q, this expression has the form of the GW ap- 
proximation [56] except for two crucial differences: (i) the effective interaction V^^,{q, ^2) in- 
cludes vertex corrections and is therefore more general than the screened interaction W{q, Q) 




'p'+qiO' p'+Qn-iCJ^ 




G,,.(p') 



FIG. 2: Structure of the diagrams included in the approximate evaluation of /po-,p'o-'(0)- The 
number of irreducible interaction blocks is n > 2, and there is an integral over the internal four- 
momenta (7i, ...,g„_i. Notice that these diagrams are reducible in the particle-hole channel that 
carries four-momentum p — p' , but irreducible in the particle- hole channel defined by the two 
external legs, which carries zero four-momentum. The irreducible interaction blocks are further 
analyzed in Fig. 3. 

between test charges that appears in the GW approximation; and (ii) the expression (|14p 
involves an undetermined integration constant which must be fixed by independent means, 
for example by requiring that T^lkp, 0) reproduces the correct value of the chemical potential 
as determined from QMC data. An analytic continuation procedure allows one to recast the 
time-ordered self-energy in Eq. ()14p into a retarded self-energy given by the sum of SX and 
CH contributions as in Eq. (|^. 

In their derivation of the "local approximation" in Eq. (J13j) Ng and Singwi sorted the 
diagrams that contribute to the irreducible electron- hole interaction into three classes: (1) 
diagrams that are reducible in the "crossed" particle-hole channel, carrying momentum q and 
energy Q; (2) diagrams that are reducible in the particle-particle channel; and (3) diagrams 
that are irreducible in the particle-particle channel and in both particle-hole channels. The 
diagrams of class (1) are shown in Fig. 2. Notice that these diagrams are expressed in terms 
of three building blocks as shown in Fig. 3, which are irreducible in the "crossed" particle- 
hole channel. In block (a) the particle and the hole have parallel spin orientations which 
are conserved as they scatter against each other. In block (b) the particle and the hole have 
parallel spin orientations, which are reversed as a result of the scattering process. Finally, 
in block (c) the particle and the hole have opposite spin orientations. 

In order to make further progress we now approximate these irreducible blocks by func- 
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(a) 



(b) 



(c) 




p+qii 





p'+qii 



FIG. 3: The three irreducible blocks out of which the "railroad diagrams" of Fig. 2 are constructed. 




p,-o 



p+q,-a 




p,a 



p'+ q, o 




P,CJ 



p+q,-a 



FIG. 4: Graphical illustration of the relationship between different components of the irreducible 
electron-hole scattering amplitude in a spin-unpolarized state. 

tions of the electron-hole momentum q = k — k'. The form of these functions is determined 
by requiring that the same approximation, when applied to the evaluation of the diagrams 
for the density- density and spin-spin response functions, yields Eqs. (jHl) and ®. For the 
two diagrams 3(a) and 3(b) this is accomplished by setting 



v^[l - a.(q)] 



(15) 



where a and r are the spin orientations of the particle and of the hole before and after 
scattering (the labelling is explained in Fig. 3). For diagram 3(c) a simple reasoning based 
on isotropy in spin space for the paramagnetic state leads to 



Tcrcr.acr jaa.aa jcrcr.crcr 



-2t;qG_(q) 



(16) 



as one can see from Fig. 4. 

The integrals over the internal four-momenta in Fig. 2 can be carried out analytically if 
the intermediate Green's functions are replaced by noninteracting Green's functions, giving 
a result proportional to [xo{p — p')]"' where xo is the Stern function. The whole series of 
diagrams can then be summed algebraically yielding 

+ 2KG_(q)]2x5(q,fi)5.,-.'- (17) 



11 



The subscripts + and — denote as usual the charge and the longitudinal spin channels. The 
density-density and spin-spin response functions are expressed in terms of the Stern function 
and of the many-body local field factors according to Eqs. (0) and ©• 

The evaluation of the remaining terms, i.e. the particle-particle ladder diagrams and the 
fully irreducible diagrams, is considerably more complex. The only simple diagram is the 
bare exchange interaction with momentum transfer q, which yields the Hartree-Fock self- 
energy. All other terms are dropped in the possibly naive hope that they are smaller than 
the terms retained in Eq. ()17p. Keeping only the bare exchange interaction and combining 
it with Eq. ()17|1 we arrive at an effective interaction of the form 



V.f{q,n) = 2[v^G4ci)fxsici,n) 



(18) 

U 

Inserting this in Eq. (|T^ and repeating standard analytical transformations [^ one easily 
recovers the expressions and (jTj) for the screened-exchange and Coulomb-hole contribu- 
tions to the self-energy. 

Let us emphasize again that the result that we have obtained by a diagrammatic method 
rests on Eq. (fT^ for the self-energy (apart from an additive constant) and on the use of the 
effective interactions shown in Eq. ()18|1 . It should be evident from the above derivation that 
no diagram for /kwo-.k'w'o-' (0) has been double-counted. Rather, many diagrams have been 
dropped, but the result for the self-energy remains to be adjusted a posteriori by fixing the 
addivite constant through the correct value of the chemical potential. 

Let us see, on the other hand, what would have happened if we had applied the local 
approximation directly to the self-energy starting from the alternative exact expression, 

where A is the proper vertex function and e is the regular dielectric function (see Fig. 5). 



Within the local approximation one finds 



la 



26 



1 
1 + t;k-k'G'+(k - kOxo(k - k', cj - uj' 



^ko^.k'o;' - ^^z. TTTT: 1,/N,, n, "TT-. ^ ' (20) 



so that this route to the self-energy includes only the contribution of charge fluctuations 
but misses completely that of spin fluctuations. The root of the difficulty obviously lies 
in the fact that the local approximation for the vertex function is not good enough to 
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W(q) 




Ga(p+q) 



FIG. 5: Diagrammatic representation of Eq. (jT^. Here 1/F(q, 17) = fq/e(q, 17) is the usual test 
charge-test charge screened interaction. 

capture the contribution of spin-density fluctuations. On the other hand, the dependence of 
/ on spin fluctuations is manifest in the terms proportional to G^ in Eq. (fTTjl . This is the 
main physical reason why it is better to apply the local approximation to the differential 
relation p2|) than to the integral relation (fT^ . In fact, all quasiparticle properties of our 
present interest depend on relative variations of the self-energy, i.e. on 5S rather than on 
the absolute value of S. 

B. Renormalized Hamiltonian approach 

In this Section we discuss the derivation of Eqs. and ((Zj) from the point of view of 
an effective renormalized Hamiltonian for the low-energy degrees of freedom of the electron 
liquid y,y,bl|. 

We start by dividing the Hilbert space of the EG Hamiltonian Heg reported in Eq. ([T} 
into a "slow" sector (5a) and a "fast" sector (J-'a), assuming the existence of the Fermi 
surface aX k = kp. S^ contains only plane- wave states with wavevector k close to the Fermi 
surface, i.e. such that \k — kF\ < A where A is an arbitrarily small cutoff. JF^ contains all 
the other states. We correspondingly introduce "slow" and "fast" creation and annihilation 
operators which operate in these two sectors, 

J Sk,„ kG^A, 
Ck,a = < . (21) 

y fk,a, k G J^A • 

Our aim is to derive an effective Hamiltonian for the slow sector which contains only the 
■Sk.cr operators by integrating out in a reasoned manner the /k,cr degrees of freedom. 
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TYeg is first rewritten using the Sk,o- and /k,(T operators, 

'^^EG = Hs + Hi + Tisf ■ (22) 

The first term is 

Hs = J2^^^ij^'- + ^ 5Z^q 5Z 5Z 4i+q,ai42-q,a2^k2,a2Ski,ai , (23) 

k,o- qT^O ki.cri k2,cr2 

where all wavevectors belong to 5a. Note that the second term in Tig, which represents the 
direct interaction between slow particles, tends to zero faster than the kinetic energy term 
in the limit A — i> 0. Thus this term will be treated by first-order perturbation theory below. 
The second term in Eq. ()22j) is the Hamiltonian for the fast sector. 



^f - 2^ ^k/k,<7/k,a + ^ Z^^q Z^ Z^ /ki+q,ai/k2-q,a2/k2,<T2/ki,(Ti (24) 

k,o- qj^O ki,(Ti k2,cr2 

where all wavevectors belong to JF^. 7i{ tends to the full EG Hamiltonian in the limit A ^ 
and this property will be useful in what follows. 

The third term Tigf describes the interaction between the slow and the fast particles: this 
term is the sum of fourteen different terms, but here we assume that the relevant operators 
are those which separately conserve the number of particles in the two sectors, i.e. terms 
of the type ■5|^^_|_ ^^s|^^ ,^Sk2,o-2/ki,o-i (trilinear in the field operators of either slow or fast 
particles) or sl.^j^ ^^sl.^_^^^fk2,a2f\ii,ai will be dropped. Four terms are left within this 
assumption. 



"sf — 2S ^-~^^'^ ^-~^ / . ■^ki+q.o-i/k2-q.o-9/k2,o-2'^ki,o-i 
q^O kijcri k2,o-2 



"*" 23^^"^^ ■' ■^ki+q.o-i'^k2-q.o-^'^k2,cr2/ki,cri 

qT^O kijcri k2,o-2 

"*" 2S 2-^^"^ 2-^ / . ■^ki+q.(n/k2-q.(T^'^k2,CTiJki,CT2 
q^O kijcri k2,o-2 

+ 2^ 2-^^^ 2-^ Z^ /ki+q,(Ti'^k2~q,o-2/k2,o-iSki,cr2 (25) 

q^O kijcri k2,o-2 

where ki and ki + q belong to S\, and k2 and k2 — q belong to JF^. By simple algebraic 
manipulations Tigf can be written as 

"sf = ~Q 2-^ 2-^ / J ^q'^ki+q,o-i/k^~q.o-2.^k2,o-2-^ki,o-i 
qT^O ki,(Ti k2,(T2 

~ ~d Z-^ 2-^ Z^ ^'^2-ki-q'5ki+q,o-i/k2-q,cr2/k2,cri'5ki,(T2 • (26) 



q^O ki,(Ti k2,o-2 
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The first (direct) term in Tisf describes the Coulomb interaction between slow and fast 
particles and can then be expressed in terms of density fluctuations in the two sets of 
particles. The second (exchange) term describes an exchange process in which a slow particle 
replaces a fast particle and vice versa. Note that the two particles can have opposite spins, 
i.e. 0"! = — o"2, and when this is the case a net spin angular momentum is exchanged between 
slow and fast particles. This in practice means that any attempt to write Tigf in terms of 
collective variables must also involve an interaction between slow and fast particles mediated 
by spin fluctuations. 

All these arguments bring us to a second crucial approximation: we treat ?isf in an average 
sense by writing these microscopic processes in terms of interactions between density and 
spin-density fluctuations in the two sets of slow and fast particles, 

q^O k.o-i 

+ ^ E ^s(q) E S-. ■ (4-q,.. W.,.. Sk,.,) , (27) 

where n^ = Ek,<7 /k-q.a/k,^ and Sq = Ek,<7,<T' /k-q,<7 Haa' hy are, respectively, the density 
and spin-density operators for the fast sector. The effective interaction potentials fc(q) 
and fs(q) must include both the exchange and the correlations effects. This requirement 
can be fullfilled in an approximate way by means of local field factors G+(q) and G-(q) 
by taking fc(q) = ^^[1 — G'+(q)] and fs(q) = —VciG-{q). Note that we are again using 
frequency-independent local field factors for the reasons already indicated above. 

We now carry out a unitary transformation that eliminates the interaction Tigf be- 
tween slow and fast particles to second order in its strength 25|. We search for a 
Hermitian operator Ga which maps the original Hamiltonian into a new Hamiltonian 
H' = exp (iG)A)'HEGexp (— zGa) having the same eigenvalues but transformed eigenfunc- 
tions. The generator Oa is at least of first order in the strength of the interaction between 
the two sectors. The transformed Hamiltonian can then be expanded according to 

n' = 7T;EG + ^[0A,^EG]-^[eA,[eA,^EG]] + - 

where we have dropped the commutator [6a, [G)A,7Ysf]] because it is of at least third order. 
The interaction term Hsf is eliminated by choosing Ga as the solution of the operatorial 
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equation 



z[eA,7^s + Hf 



-n 



sf 



(29) 



Finally, by averaging over the ground-state |0) of the Hamiltonian Tif at energy Eq, we 
obtain the effective Hamiltonian for the low-energy degrees of freedom of the electron liquid, 



T^QP = (o|7^'|o) = E, + n. + _(o|[e^,7^,f]|o) . 



(30) 



Obviously Eq does not play any physical role and will be dropped from now on. 

The operatorial equation (j^^ can be solved for Ba once the commutator of Ba with the 
interaction term in Ti^ is dropped, the justification being that this commutator vanishes upon 
averaging on the ground-state of the fast sector. A lengthy but straightforward calculation 
yields 



q^O 



ki,o-i k2,(T2 



^Xl^s(q) Y^ Y^ ^s(kl,k2,q)4i+q,<xiW-l-25ki 



ci 



q^O 



E 



ki,o-i k2,(T2 



'^k2-q,Ti ['''J-riT2'^k2,r2 



(31) 



where we have defined 

Mc(ki,k2,q) = ^5^ 

n 

and 

Ms(ki,k2,q) = -5^ 



(n|riq|0)p 



[n 



n,|0)p 



£^k2 



£k 



2-q 



Eq — En ^ki — ^ki+q + En — Eq 



l(^l^.,q|0)P 



(^l'^.,q|0)P 



£^k2 



£k 



2-q 



Eq — En £:ki — ^ki+q + En — Eq 



(32) 



(33) 



Here \n) is an exact excited eigenstate of the fast-sector Hamiltonian with eigenvalue En- 

A key point is that Tiqp as written in Eq. ()3H1 is not normal-ordered with respect to 
the vacuum, i.e. contains a self-interaction term that needs to be subtracted. The physical 
meaning of this term is clear: a slow particle moving through the "medium" of fast-moving 
particles creates a local polarization which in turn acts back on it. To subtract this self- 
interaction term we just proceed by normal ordering Tiqp with respect to the vacuum, with 
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the result 



^QP = Yl 



k.o- 



^k 



1 ^ r'^ duo' 



q 



S "^^ Jo vr Ak,q - to' 



sIJk< 



+ ^5Zt^'l + ^c(q)Xc(q,Ak,q)] XI 5Z C+q-lC-q.aa^^ka.a^Ski.ai 
q^O ki,(Ti k2,cr2 

+ ^5Z^s(q)Xs(q, Ak,q) XI Zl 5Z W-l-2 ■ Wnra^i+q.ai^a-q.n^ka.r^Ski, 



QT^O kijCTl k2,(T2 Tl,T-2 

(34) 

where Ak,q = (e^k — £^k-q)/^- The sums over the eigenstates of Hi in Eqs. flH^ and (jH^ have 
been carried out using the identity 

S ^ Ek- £k-q - [E„. - Eq) tt Jo (^k - £k-q)/^ - t^' 

and a similar identity for S'^^q. Here Xc(q; i^) and Xs(q! ^) are the response functions of the 
"medium" , which asymptotically tend to those of the full EG for A ^ 0. 

The QP Hamiltonian Tiqp as written in Eq. (jMj) has a very clear physical meaning. 
It is the Hamiltonian for a gas of weakly interacting slow-moving particles with a single- 
particle dispersion relation shifted by a self-interaction term (the so-called Coulomb-hole 
shift) generated by the normal ordering with respect to the vacuum. The reason for being 
weakly interacting is that the quartic term in Tiqp tends to zero faster than the kinetic- 
energy term in the limit A -^ 0. We can thus calculate the quasiparticle energy within 
first-order perturbation theory, with the result 

^QP (k) = £k + ^sx (k) + ^cH (k) . (36) 

Here the first term is the bare single-particle energy and the second term, which has been 
generated from the normal ordering with respect to the Fermi sea, is formally an ordi- 
nary exchange self-energy calculated with a dynamically screened effective interaction (the 
Kukkonen-Overhauser QP-QP interaction): 

^sx(k) = - j -^ K + t;2 (q)xc(q, Ak,q) + 3i;i(q)xs(q, Ak,q)] Q{-i^.-Jh) . (37) 

The last term in Eq. (|HH|) is the Coulomb-shift generated by the normal ordering with respect 
to the vacuum, 

^ (A\ f ^^^ /"^"^c^o;' f^(q)Smxc(q,^') + 3t^s(q)^"^Xs(q,^') .oqn 

J (27r)^ Jq tx Ak,q - uo' 
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The sum of £^sx(k) and ScuiX) coincides with Eq. (jH)) calculated at the single-particle fre- 
quency u = C,k/h. 

III. LOCAL FIELD FACTORS 

As is clear from Eqs. Q-©, the local-field factors are fundamental quantities for an 
evaluation of quasiparticle properties. In this Section we introduce the static values of these 
functions, that we have chosen to calculate the real and imaginary parts of the QP self-energy 
reported in Eqs. (jH)) and Q. 

Analytical expressions are availabl e |3ll| for G+(q) and G_(q), which reproduce the most 
recent Diffusion Monte Carlo data |7|, |l^ and, as we are going to summarize below, embody 
the exact asymptotic behaviors at both small and large wave number q. Specifically, in 
the long wavelength limit our choice satisfies the compressibility and spin-susceptibility sum 
rules, 

hm G±(q) = A±-^ (39) 

with A^ = (1 — kq / k) / [r s\/2) and A^ = (1 — xp I Xs) / {'^ s^/^) ■ Here kq is the compressibil- 
ity of the noninteracting gas, n and xs ^'^^ the compressibility and the spin susceptibility 
of the interacting system, and xp is the Pauli spin susceptibility. By making use of the 
thermodynamic definitions of k and xs we can write 



(9r? r, dr^ 



(40) 



and 

(41) 



XP _,_V2 rl d^e,{r,X) 

i f s ' 



xs vr - 2 de 



C=o 



where ..(r,.C).s the corre 
of spin-polarization |10lll5| 



59| 



ation energy per particle as a function of r^ and of the degree C, 
At large g, on the other hand, the local fields of Ref.|3l| satisfy the asymptotic behavior 58, 



G'±(q)^C-^ + S±. (42) 

kp 

Here C is determined by the difference in kinetic energy between the interacting and the 
ideal Fermi gas, 
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while -B+ = 1 — g{0) and -B_ = g{0), with g{0) being the value of the pair distribution 
function at the origin 60|. 

IV. NUMERICAL RESULTS 

We turn to a presentation of our main numerical results. In Sect. IIVAI we present 
some illustrative results for the QP excitation energy and lifetime, and in Section IIV Bl we 
give our results for the QP effective mass and renormalization constant. In all figures the 
labels "RPA" , "6*+" and "G+feG^" refer to three possible choices for the local-field factors: 
"RPA" refers to the case in which local-field factors are not included, "(?+" to the case in 
which the antisymmetric spin-spin local field is set to zero {i.e. spin-density fluctuations 
are not allowed), and finally "G+feG-" refers to the full theory including both charge- and 
spin-density fluctuations. 

A. Quasiparticle self-energy 

We have computed the real and imaginary parts of the QP self-energy using Eqs. (jHI) and 
(jH)). In Figs. 6 and 7 we show the real part of the SX and CH contributions as from Eqs. (0) 
and ((Zj), evaluated at the single-particle frequency u = C,k/h and measured from their value 
at fc = kp. Note the presence of a strong dip in the CH term at a value of k {kp, say) 
which depends on r^ and on the functional form of the charge-charge local field factor. This 
is the plasmon dip, which is also present in 3D and originates from the fact that at each 
Tg there is a sufficiently high value of k for decay of an electron-hole pair into a plasmon 
with conservation of momentum and energy. Mathematically this dip arises for the reasons 
explained in Appendix A. 

In Fig. 8 we show 3f?eS|!^^.(k, u) as from Eq. (jHJ, evaluated at u; = ^k/h. There is substantial 
cancellation between the SX and CH contributions for k < kp, so that in this range the QP 
self-energy is essentially very weakly momentum-dependent. Such a function has a Fourier 
transform which is to a good extent local in real space, and this result can be viewed as 
a microscopic justification of the local-density approximation to the exchange-correlation 
potential of density-functional theory. 

In Fig. 9 we show the absolute value of Q^mSrotlk, i^) as from Eq. Q, evaluated at 
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uo = C,]i/h. This function takes a finite jump at the wave number of the plasmon dip. The 



El: 



discontinuity is pecuhar to 2D |6JJ: it is absent in 3D and arises from the fact that the 
oscillator strength of the plasmon pole is non-zero at kp (see Appendix A). The qualitative 
difference in the shape of the imaginary part of the self-energy below and above kp reflects 
the opening of a new decay channel for an electron-hole pair. 

B. Many-body effective mass enhancement 

Once the QP excitation energy is known, the effective mass m* can be calculated by 
means of the relationship 

(44) 



1 1 dS£Qp{k) 



m* h^kp dk f,^j_. 
In Sect. |n]we remarked that the QP excitation energy may be calculated either by solving 
self-consistently the Dyson equation (fTn|) or by using the OSA in Eq. (fTTj) . In what follows 
the identity 

— dfc5KeS,,t(/c,^)|^^^(^,^ + dJRe\^,[k,u)\^^^^^^ —^ (45) 

will be used, uj{k) being an arbitrary function of k. 

Using Eqs. (j44p and (j45|) with uj{k) = 6SQp{k)/h we find that the effective mass ttt.^ 
calculated within the Dyson scheme is given by 

'j}ii = ?:! (46) 

The renormalization constant Z that measures the discontinuity of the momentum distri- 
bution at k = kp is given by 

The normal Fermi-liquid assumption, < Z < 1, implies c)^3fieE^^(/c, co')|^_^ ^ < 0. 

On the other hand, using Eqs. (jl^ and (jl^ with u{k) = C,]i/h we find that the effective 
mass ^TiosA within the OSA is given by 



OSA ^ ^ /^g-j 

Of course, Eq. PHj) is a good approximant to the Dyson effective mass only in the weak- 
coupling limit as can be seen by expanding Eq. pUj) near the point r^ = 0. The normal 
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Fermi-liquid assumption with the accompanying inequahty otherwise imphes that a zero of 
the denominator in Eq. (jlSj) might occur depending on the inputs which are used to calculate 
the QP self-energy. In fact it is clear from Eq. (|46|) that a singular behavior of the effective 
mass can only occur if the renormalization constant goes to zero, i.e. the normal Fermi- 
liquid assumption breaks down. In this case the singular behavior could be interpreted as a 
quantum phase transition of the 2D EG to a non- Fermi-liquid state. 

In Fig. 10 we show our numerical results for m^ and w^qsa- The effective mass enhance- 
ment is substantially smaller in the Dyson-equation calculation than in the OSA, the reason 
being that a large cancellation occurs between numerator and denominator in Eq. ()4(ij) . 
In both calculations the combined effect of charge and spin fluctuations is to enhance the 
effective mass over the RPA result, whereas the opposite effect is found if only charge fluc- 
tuations are included - a manifestly incorrect result that neglects the spinorial nature of 
the electron. For completeness we have also included in Figure 10 the variational QMC 
results of Kwon et al. js|. The reader should bear in mind that the effective mass is not 
a ground-state property and thus its evaluation by the QMC technique is quite delicate, 
as it involves the construction of excited states. There clearly is quantitative disagreement 
between our "best" theoretical results (the "G+feG-ZD" predictions) and the QMC data. 

In Fig. 11 we show the behavior of the two terms in the denominator of Eq. ()48|) as 
functions of r^. This figure clearly shows how a divergence can arise in wiQg^: for instance, 
within the RPA the denominator in Eq. ()48|1 has a zero at r^ ~ 15.5 (see the inset in 
Fig. 11). Our numerical evidence, within all the three theories that we have studied, is that 
indeed (i) d^^j'^eYi^^^ikp^Q) is negative as it should for a normal Fermi liquid, and monoton- 
ically increasing in absolute value as a function of r^; and (ii) dk'SleTi^^^^ikF^Q) is positive 
and monotonically increasing too. Within the theory outlined in this work, which uses as 
a key ingredient the Kukkonen-Overhauser effective screening function in Eq. (j3)), the ef- 
fect of a charge-only local field is to shift this divergence to higher values of r^, while the 
opposite occurs upon including both charge and spin fluctuations. For instance, within the 
"G+feG-ZOSA" theory the divergence occurs near r^ = 5. Within the local approximation 
of Eq. (j^Uj) the situation is different 6^ and the effect of a charge-only local field is to shift 
the divergence to lower values of r^. 

In Fig. 12 we show our numerical results for the renormalization constant Z in comparison 



with the QMC data of Ref. |63|. The theory underestimates the value of Z over the whole 

21 



range of densities explored. Notice that short-range charge- density fluctuations tend to 
stabilize the normal Fermi liquid, while the simultaneous inclusion of charge- and spin- 
density fluctuations works in the opposite way. 

For the sake of completeness we have collected in Table I a summary of our numerical 
results for m* and Z at a few values of r^. In the next Section we will discuss our results 
for the QP effective mass in the light of recent experimental results and draw our main 
conclusions. 

V. COMPARISON WITH EXPERIMENTAL RESULTS AND CONCLUSIONS 



A full analysis of the published data for the effective mass of carriers in Si-MOSFET's 42 , 



^ 
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43l | would require a more complete theoretical study, mainly to account for the two-val 
nature of the material. We will focus here instead on the experimental results of Ref. 
as kindly provided to us by Dr. Zhu prior to publication. At present the data refer to the 
range 2 < r^ < 6, so that we cannot judge the performance of the theory in the weak- 
coupling regime. A quantitative comparison between theory and experiment would also 
require a refined treatment of a series of effects such as those due to disorder and to finite 
temperature. We restrict our analysis to the effect of finite sample thickness, by discussing 
how a softened Coulomb potential modifies m* against the strictly 2D results discussed in 
Sect. IIVBI and shown in Fig. 10. The expectation is that the QP effective mass will be 
noticeably smaller when a softened Coulomb interaction is at work. 

We have thus recalculated m* after renormalizing the bare Coulomb potential by means 



le GaAs/AlGaAs 
The appropriate 



of a form factor to take into account the finite width of the EG in t 
heterojunction-insulated gate field-effect transistor used in Refs. I46landl51 
renormalized potential is given by V^ = Vc^F{qd)/R, where 

with d = [h? Ksc/ (4:871711^6^12*)]^'^ representing an effective width of the 2D EG J22|- Here 
Kins = 10.9 and Ksc = 12.9 are the dielectric constants of the insulator and of the space 
charge layer, k is their average, m is the band mass in the confinement direction, and 
n* = ndejji+ 11^20/32, the depletion layer charge density n^epi being zero in the experiments 



of Ref. 



5ll The results that we obtain with the softened potential are shown in Fig. 13. A 
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caveat to keep in mind is that we have used the same local-field factors as for a zero-thickness 
2D EG in the lack of a better choice. Thus the results labelled by "(?+" and "G+feG-" in 
Fig. 13 contain the effect of finite thickness only through the renormalization of the Coulomb 
potential. We believe that the explicit dependence of the local fields on the finite width of 
the 2D EG should not change the results of Fig. 13 in a substantial manner. 

Comparing the results of Fig. 13 with those in Ref. |5l| we can draw the following conclu- 
sions: (i) the "(j+" results, at both the OSA and the Dyson-equation level, do not have the 
proper functional shape to account for the experimental data; (ii) the RPA and "G+fcC-" 
results are rather similar; and (iii) the "(?+&(?„" results, which treat charge and spin fiuc- 
tuations on the same footing, show the best performance against the experimental data. In 
fact, without the use of any fitting parameters, the "G+&G_" results within the OSA com- 
pare in a very reasonable manner with the data. The Dyson-equation results, show instead 
a relatively small and slowly increasing mass enhancement over the whole range of densities, 
as discussed for the strictly 2D case in Sect. IIVBI 

In summary, we have revisited the problem of the microscopic calculation of the quasipar- 
ticle self-energy and many-body effective mass enhancement in a 2D EG. We have performed 
a systematic study based on the many-body local-fields theory, taking advantage of the re- 
sults of the most recent Diffusion Monte Carlo calculations of the static charge and spin 
response of the EG expressed through static local-field factors. We have carried out exten- 
sive calculations of both the real and the imaginary part of the quasiparticle self-energy. 
We have also presented results for the effective mass enhancement and for the renormal- 
ization constant over a wide range of coupling strength. In this respect we have critically 
examined the merits of the on-shell approximation versus the Dyson-equation calculation. 
Depending on the local-field factors, the OSA predicts a divergence of the effective mass 
at strong coupling and a solution of the Dyson equation is necessary in order to obtain a 
well behaved effective mass. The comparison with the experimental data of Ref. |^ shows 
that the simultaneous inclusion of charge- and spin- density fluctuations beyond the Random 
Phase Approximation is crucial in accounting for exchange and short-range correlations, and 
can lead to substantial corrections at low carrier densities. A possible role of dynamic cor- 
relations, entering through the frequency dependence of the local-fleld factors, remains to 
be examined. 
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Appendix A. Details on the explicit calculation of the real part of the CH contri- 
bution 



Using Eq. (fTj) we find that the real part of the CH term evaluated at u = ^k/h is given 



by 



^eEcnik, uj) 



l'^=Ck/?i 



VaV 



(2vr) 



2 ^q' 



dQ '^m[e-\q,Q)] 







TT 



^k/^ - ^k+q/^ - ^ 



(50) 



The angular integration can be performed analytically, with the result 

Jo ^k/h-^k+Jh-n y/[Q + %V(2m)]2 - h^k^q^w? 

where r2jnin(^) = —hq^/{2m) + hkq/m. In carrying out the frequency integration care must 
be taken to include the contribution from the plasmon pole Qp\. Using the expression for 
the imaginary part of the charge-charge susceptibility near fipi, 

-1 



3mxc(q,^) = 7rwc^(q) 3f^exo(q,fi)| 



f7=n„ 



93f?exo(q,fi) 



dn 



n=n„ 



6{Q - fipi) 



(52) 



we find that the real part of the CH term is given by 



afJeScHCk, u;) 



qdq 






2n ^[f]pi + V/(2m)]2 - H^k^q^ 



95Rexo(q,fi) 



dn 



-\ -1 



n=f7„ 



6(r2pi - f2mm(fc)) 






P 



^ 'up 



dfi i;^(q)Smxc(q, ^) + 3t;s(q) S"^xs(q, fi) 



max[0,r2min(fc),f^low] ^ 



v/[fi + V/(2m)]2 - ;i2p^2/ 



r/i^ 



(53) 



hq^ /{2m) ± fiqkp/m are the 



Here gc marks the onset of Landau damping and i^up(iow) 
upper and lower edges of the electron-hole continuum. 

The range of the momentum integration deserves special attention. In the first term in 
Eq. ()5H|1 ■ due to the step function the range of g- integration is determined by the intersections 
between r2min(fc) and Vt^y (see Fig. 14). In Sect. llVXl we have introduced the r^-dependent 
wave number kp. this is the wave number k at which fijnin(^) is tangent to Vt^i. There are 
two cases: (i) for k < kp there are no intersections and thus the range of g-integration goes 
from to qc', and (ii) for k > kp there can be either one (gi) or two intersections (gi.2), so 
that the range of integration is [0, gi] or [0, gi] U [g2, qd, respectively. It is the crossover from 
condition (i) to (ii) that leads to the plasmon dip in 9fJeScH(k, i^)L=f^/?i- 
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In the second term in Eq. ()53p the range of g-integration runs up to g = +00 and this gives 
rise to the logarithmic divergence mentioned in Ref. [591 As aheady discussed in Sect. Ull 
what matters are self-energy differences, which are free of singularities. Numerically we deal 
only with the finite quantity ^eT,cR(k,uj)\^^^ ,^ — SchI^f^O). 



Appendix B. "Line+residue" decomposition 



il 



In this Appendix we discuss a mathematical 
energy, first introduced by Quinn and Ferrell 
literature 6J|. This amounts to writing 



Si.et(k,u;) = Sx(k) + Sii„e(k,u;) + Sres(k,u;) 
Here the first term is the Hartree-Fock self-energy 



equivalent decomposition of the QP self- 



17j . which has been often employed in the 



(54) 



Sx(k) 



(55) 



-2e2fci.E(P)/7r (^ < 1) 

-2e^kF k [E(l/P) - (1 - 1/P) K(l/P)]/7r {k > 1) 

where k = k/kp and K(a;), E(x) are complete elliptical integrals of the first and second kind, 
respectively. The second term in Eq. ()54|). which is purely real, is given by 



Eiine(k,t^) 






,2"q 



2^ 



e(q,ifi) 



UJ + iVL- ^k+q/^ ' 



(56) 



Finally, the third term is the so-called "residue" contribution. 



Sres(k,Cj) 



d\ 



[e(a;-a+q/^)-0(-ek+q/^)]- (57) 



(27r)2 ''[£(q,a;-aV^) 
Within this decomposition it is the "line" contribution which needs to be regularized for a 
ultraviolet divergence. 

As a check of our numerical results obtained by means of the SX-CH decompositions, 
we have recalculated the QP self-energy, effective mass, and renormalization constant by 
this alternative route. This turned out to require a substantially harder numerical effort. 
For completeness we summarize in Figs. 15 and 16 our results for the "line" and "residue" 
terms. 
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TABLE I: QP effective mass and renormalization constant of a paramagnetic 2D EG. The arrows 

refer to the divergence in the OS A as explained in the main text. 

Ts Various calculations m^/m mQgj^/m Z 

1 RPA 1.022 1.033 0.670 
G+ 0.972 0.961 0.710 
G+kG- 1.026 1.040 0.658 

2 RPA 1.082 1.168 0.526 
G+ 1.004 1.007 0.585 
G+kG- 1.144 1.349 0.486 

3 RPA 1.121 1.322 0.444 
G+ 1.030 1.061 0.510 
G+&iG- 1.247 2.026 0.391 

5 RPA 1.167 1.696 0.349 

G+ 1.066 1.172 0.419 

G+&cG- 1.410 /■ 0.289 

10 RPA 1.215 3.650 0.244 

G+ 1.100 1.415 0.311 

G+&cG^ 1.834 /■ 0.166 
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FIG. 6: The real part of the SX contribution to the retarded self-energy (in units of Ryd) evaluated 
at a; = Ck/K as a function of k/kp for r^ = 1 (top panel) and 5 (bottom panel). 
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FIG. 7: The real part of the CH contribution to the retarded self-energy (in units of Ryd) evaluated 
at a; = Ck/K as a function of k/kp for r^ = 1 (top panel) and 5 (bottom panel). 
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FIG. 8: The real part of the retarded self-energy (in units of Ryd) evaluated at u; = iia/fi, as a 
function of k/kp for r^ = 1 and 5. /iQMC is the chemical potential from the QMC ground-state 
energy. 
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FIG. 9: The absolute value of the imaginary part of the retarded self-energy (in units of Ryd) 
evaluated at u; = ^k/^, as a function of k/kp for r^ = 1 (top panel) and 5 (bottom panel). 
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FIG. 10: Effective mass enhancement as a function of r^ for < r^ < 10. The inset shows an 
enlargement of the results for r^ < 1. The lines show the results from Eq. (|4(>] 1. while the symbols 
(except for the dots) are from Eq. (|48|) . The QMC data (dots) are from Ref. |a. 
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The inset shows this divergence occurring within the RPA at r^ ~ 15.5. 
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FIG. 12: Renormalization constant Z as a function of r^ for < r^ < 10. The QMC data have 
been taken from Ref. |23- 
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FIG. 14: Illustrating the integration range relevant to Eq. (|53|) . with fi in units of Epjh. The 
shaded area represents the electron- hole continuum. ^xa\\\{(i) is shown for three values of k (dashed 
curves): from top to bottom, k = 2>.{)kF, kp, and 0.8kp with kp ~ 2.32kF. The plasmon dispersion 
relation ^pi{q) (dotted curve) is also shown. 
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FIG. 15: The exchange plus regularized "line" contribution to the retarded self-energy (in units of 
Ryd) evaluated at w = ^k/^, as a function of /c/fci? for r^ = 1 and 5. 
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FIG. 16: The real part of the "residue" contribution to the retarded self-energy (in units of Ryd) 
evaluated at u; = ^k/^, as a function of k/kp for r^ = 1 (top panel) and 5 (bottom panel). 
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